o 
o 

(N 



cd 



*T3 

O 
o 



X 
S3 



Adsorption of a semiflexible polymer onto interfaces and surfaces 



Semjon Stepanow 

Universitat Halle, Fachbereich Physik, D- 06099 Halle, Germany 

(February 1, 2008) 



2 ' Abstract 

We consider the adsorption of a semiflexible polymer chain onto interfaces 
and surfaces by using the differential equation of the distribution function 
G(R, L) of the end-to-end distance R, which is associated with the moment 
c/2 ' expansion of the latter. We present the results of the approximative treatment 



consisting of taking into account the 2nd and 4th moments in the differential 
equation for G(R,L). The essential features of adsorption of the semiflexible 
polymer are: i) the existence of a new local length scale, which results in 
two-exponential decay of the monomer density of adsorbed polymer; ii) the 
binding of the semiflexible polymer is weaker than that for flexible one for 
both interface and wall. The approximative theory presented is restricted to 
the regime of weak adsorption, where the effect of the rodlike behavior of the 
polymer on small scales is weak. 
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I. INTRODUCTION 

Adsorption of polymers onto surfaces and interfaces plays a major role in polymer adhe- 
sion, stabilization of colloidal suspensions, and also in processes including biological macro- 
molecules such as protein adsorption on membranes. For many applications (including 
biopolymers) it is important to consider the effect of the stiffness on the adsorption. Over 
many years there is current interest on adsorption of semiflexible polymers in the literature 
|jlj]- ||TD| (and citations therein). Besides the interest in problem of adsorption of semiflexible 



polymers there is also large interest in bulk properties of semiflexible polymers |TT[- [13 
Models of semiflexible polymers have also applications in different topics |)14J| - |I5| ]. 

In the present article we present a theory of adsorption of a semiflexible polymer general- 
izing the approach based on the analogy between statistics of polymer chains and trajectories 
of a Quantum Mechanical particle at imaginary time |16[]. According to this analogy the 
adsorption of a flexible Gaussian polymer can be mapped onto the problem of binding states 
of a quantum mechanical particle in an external potential well ]T7|]. The key point of our 
approach consists of taking into account the higher moments of the end-to-end distribution 
function of the semiflexible polymer G(R, L) to derive a differential equation for the latter, 
which generalize the Schrodinger type equation of the flexible polymer. In the case of the 
piecewise constant potential the adsorption problem can be mapped on a quantum mechan- 
ical one: the solutions of the differential equation for G(R, L) in regions, where the potential 
is constant, have to obey the boundary conditions at jump points of the potential and in 
infinity. The binding energy and the wave function of the ground state are obtained as a 
result. 

The article is organized as follows. Section II introduces to the formalism. Section III 
and IV consider adsorption onto a symmetric interface and a surface, respectively. Section 
V contains our conclusions. 



II. FORMALISM 

The distribution function of the end-to-end polymer distance of the Kratky-Porod model 
18] is given by path integral as follows 



G(R - Ro, L) = N J Dr(s) J exp(jk(R - R„ - £ ds^-) J[ Mi^-f - 1) x 

where l p is the persistence length. The product over s in Eq. (P takes into account the local 
inextensibility of the polymer chain. If the polymer interacts with an external field, the 
potential energy — / V(r(s)) should be added in the 2nd exponential of Eq. ([[]). Instead of 
G(R, L) it is convenient to consider its Fourier transform G(k, L) 

G(k,L)= f d 3 Rexp(-ik(R-K ))G(R-R ,L), (2) 



which we represent in terms of the moments [] as follows 

i-L °° 



G(k,L)=exp(- f dsJ2^n(s)(k 2 ) n ) (3) 

J o „=1 

Eq.(^|) enables one to derive the following differential equation for G(R, L) 

9 °f L L) - E(-l) n+ W(i)A"G = 0, (4) 



n=l 



where A = V 2 is the Laplace operator. Eq.([|) is exact for a free polymer chain. In the 
presence of interaction we add the term V(R)Cr(R, L) on the l.h.s. of Eq.(|4]). In contrast 
to the free polymer we cannot now affirm that Eq.(f|) is still exact. The reason is that the 
distribution function of the end-to-end distance G(R, L) does not possess the Markovian 
property in R space. However, it is easy to see that in the formal solution of Eq.(^) the 
potential V(R) appears correctly via a Boltzmann factor. It also follows from Eq.@ that as 
it should be the higher moments /J,2n{L) suppress the number of conformations contributing 
to G(R,L). Thus, Eq.(f|) (see also @) is expected to take into account in a correct way 
the statistical properties of the semiflexible polymer and its interaction with the external 
potential. 

The 2nd and the 4th moments of the distribution function of the end-to-end polymer 
distance of the free polymer chain were computed for the Kratky-Porod model as ||18||- pifl 





Ri = 2/ 2 (exp(-^) + j - 1), 
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RA = 8l ^ -^f + yf - 4ex p(-r) - r exp( ~r } + Yi ex p(- 3 r))- ( 6 ) 

vJfc^j \J Lp — . / L'p Lp Lp — / Lp 

These moments are related for large L according to i? 4 = 5/3 (-R 2 ) 2 (for space dimension 
d — 3), while for small L the relation reads i? 4 = (R 2 ) 2 - Demanding that R? and i? 4 given 
by Eqs.(§-|6|) are identical with those computed from Eq.(|4]) we find ^(L) and ^a(L) as 

(7) 
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H A (L) = ^(ll/ p + 3/ p exp(-f)-24Lexp(-f) + / p exp(-3-)-15/ p exp(-2-)). (8) 

ip ip 

Neglecting all moments in Eq.([|) besides the 2nd one gives the theory of a flexible polymer 
l6fl . Notice that Eq.(|) with ^(L) given by Eq.flT]) yields the Gaussian distribution function 



of the end-to-end distance with the correct 2nd moment given by Eq.(|5|). However, Eq.([7p 



1 Exactly speaking H2n{L) are cumulants of the moments. 



with the exact fi 2 (L) does not give anything new on adsorption, because the adsorption 
problem is related to the large L behavior of Eq.(|7|). The terms in Eq.(|4]) associated with 
derivatives A 2 G, A 3 G, ... are responsible for the deviation of G(R,L) from Gaussian 
distribution and are due to the bending energy and the local inextensibility of the polymer. 
Notice that despite the presence of the high-order terms in (^) the distribution function 
G(R, L) associated with Eq.@ tends for large L to a Gaussian function. The use of Eq.(f|) 
in computing the end-to-end distribution function G(R, L) |22|| - j£3] by taking into account a 
few moments has the deficiency that the latter becomes negative for finite L at values of R, 
R/L ~ 1.2. However, this value of R/L corresponds to the full stretching of the polymer and 
is thus outside the practical relevance. Thus, this deficiency does not immediately concern 
the adsorption problem, since in the latter only a few modes are relevant, while the negative 
part of G(R, L) is due to summation over all modes. The first surprise, which follows from 
Eq.(||), is that /^(L) does not vanish for large L, //4(oo) = ^11/135. The nonzero value of 
[14(00) originates from the preasymptotic terms in R? and R A . This property is important 
for adsorption of a semiflexible polymer. 

To study the adsorption of a semiflexible polymer based on Eq.@ we add on the left- 
hand side of Eq.(f|) the term U(z) G, where U(z) = x@{~ z ) — u 9(z)9(w — z) is the potential 
associated with an asymmetric interface, and is shown in Fig.l. Adsorption onto an impen- 
etrable surface can be obtained from U(z) in the limit of large x- Integrations in Eq.(^j) over 
the transversal coordinates gives an equation containing only the dependence on R z = z. 
The three dimensional character of the initial problem remains contained in the particular 
values of [i 2 and ^4. Thus, the investigation of adsorption of a semiflexible polymer reduces 
to the study of the differential equation 

^ - /2 2 d 2 z G + ^3\ G + U(z) G = 0. (9) 

Having in mind to study Eq.@ for large L we have replaced ^ 2 {L) and /^(L) in Eq.@ by 
their asymptotic values fi 2 (oc) and ^4(00). 

In solving Eq.(|9]) we follow the quantum mechanical method of treating the problem 
of bound states in a potential well described by a piecewise constant potential |23j|. The 



solution of Eq.(|) in the regions A (z < 0), B (0 < z < w), and C (w < z), where the 



potential U(z) is constant, is chosen in the form G(L,z) = exp(— eL) exp(±kJ fj, 2 /2n 4 z) 
where k takes two values hi and k 2 given by 



fc^Vl + v^, fc 2 = v /1 -v / ^i' (10) 

with ei = — 1 — -^f(e — it), and u being the value of the potential (y, — «o, 0) in A, B, and 
C, respectively. Henceforth we will measure the energy and distances in units of /i|/4/i4 



and i/2/i 4 /yU 2 , respectively. Notice that taking into account higher derivatives in Eq.(|D will 
results in more eigenmodes fcj(— e), which contribute to the space modulation of the wave 
function G(L, z). The wave vectors k\ and k 2 are real if e\ is negative and smaller than —1. 
We will assume in the following that the energy obeys the condition, — e < 1, which defines 
the regime, which we interpret as weak adsorption. We chose the solution in A and C as 

i/j a (L,z) = exp(-eL)(aiexp(fc ai i2;) + a 2 exp(k ay2 z) (11) 



ipc{L,z) = exp(-eL)(ciexp(-k Cy iz) + c 2 exp(-k C)2 z), (12) 

where k c = k a for a symmetric potential. Expecting that in the region B the energy obeys 
the condition, — e > Uq, we obtain that k\ is real while k 2 is imaginary, so that the solution 
in B can be written in the form 

ip B (L,z) = exp(-eL)(&iexp(/c M z) + b 2 exp(-k b:i z) + b 3 cos(k b>2 z) + b 4 sm(k b:2 z)). (13) 



As in the case of a flexible polymer [[H]] the density profile of the monomers in the approxi- 
mation of the ground state dominance is proportional to | ip(z) | 2 , while the wave function 
itself gives the distribution of the chain end, and, thus, should be positive quantity. In 
the following we will consider the adsorption of a semiflexible polymer onto a symmetric 
potential well, and onto an impenetrable surface, separately. 

III. ADSORPTION ONTO A SYMMETRIC INTERFACE 

The coefficients ai, ..., c 2 in Eqs.(|TT|-|r3|) have to be computed from the boundary con- 
ditions, which consist in continuity of the function ip{z) and its three derivatives at the 
boundaries z = and z — w. The boundary conditions yield a linear homogeneous system 
of eight equations, which have a nonzero solution, if its determinant A ac is zero. The con- 
dition A ac = gives the energy eigenvalues. In the vicinity of the localization transition 
(small Mo) the energy eigenvalue and the amplitude in Eqs. ( |TT| - |T3| ) are obtained respectively 
as 



2„,2 



W U 



V2 



-e = _JL(l-^ Wi0 + ...), (14) 

ai = - (u /8)(l - exp(- v / 2w)), a 2 = 1, b x = (u /8) exp(- v / 2w), b 2 = u /8, b 3 = 1, 
6 4 = (y/2/4)y/uo w , c\ — («o/8)(l — exp(\/2w)), c 2 = 1. The prefactor in front of the wave 
function is determined from the normalization of the latter. The leading term in Eq.(|HD 
coincides exactly with the energy obtained for a flexible polymer after using the same units 
(the prefactor 1/8 will become 1/4, if lengths are measured in units of fi 2 ). It follows from 
Eq. ( |14"D that the energy level for the semiflexible polymer is higher in comparison to that of 
the flexible polymer, so that the binding of the semiflexible polymer is weaker in comparison 
to the flexible one. This result can be interpreted as follows. Let us consider a small piece of 
the polymer with number of segments n, which is in contact with the interface. The size of 
that piece is proportional to n for a semiflexible polymer and to y/n for a flexible one. Thus, 
the number of contacts of the flexible polymer with the interface is larger, which results 
in a larger energy win. The consequence of the smaller binding energy of the semiflexible 
polymer is that the space distribution of the semiflexible polymer is broader in comparison 
to that of the flexible polymer. This finding is in agreement with the results of van Eijk and 
Leermakers 0. 

However, the most striking difference in comparison to adsorption of a flexible polymer 
is the appearance of the 2nd length scale £ s f ~ 1/fci, which at the localization transition is 
of order of magnitude of the persistence length l p . This length is due to the fourth derivative 
in Eq.(|4]). It is finite at the localization transition, where the localization length £ ~ l/k 2 is 
large. 



Notice that the amplitudes b\, b 2 , 01, and c\ vanish at the localization transition. The 
local length £ s f increases with increase of the depth of the potential well Uq and converges 
towards the localization length £, which decreases with increase of w . The increase of £ s f 
means that the fraction of monomers associated with local order increases. The both lengths 
C, s f and £ approach each other, if e approaches the value —1. In increasing the depth of the 
potential well, the energy tends to the value —1, while the wavenumbers k\ and k 2 approach 
the value 1. It appears that at — e = 1 the wave function becomes exactly zero. The 
threshold value u crit corresponding to — e = 1 is obtained for w = 1 as u sa = 4.24. Such a 
state does not exist, so that the description of adsorption by taking into account the fourth 
moment in Eq.(|9|) breaks down at — e > 1. For u > u sa the rodlike character of the polymer 
chain is expected to play a dominant role. The higher moments in Eq.(||) have to be taken 
into account to describe the adsorption for u > u sa . 

IV. ADSORPTION ONTO A SURFACE 

The piecewise constant potential, which models the adsorption onto a wall, is obtained 
from that shown in Fig.l in the limit x ~ > °°- To ensure that the wave function is zero at 
z = we take the latter in the region B as 

ijj B (L,z) = exp(-eL)(&isinh(/c 6)1 ,2) + b 2 sm(k bj2 z)). (15) 

The wave function in the region C is chosen according to Eq . (fL2"l) . The coefficients b\, b 2 , and 
C\ are expressed through the coefficient c 2 by using three boundary conditions at z = w. c 2 
can be fixed due to normalization of ip(z). The fourth boundary condition gives the energy 
eigenvalue, which is restricted to the condition e < — 1. The threshold value for the depth of 
the potential well is obtained for w = 1 as u t = 7.8. The amplitudes a%, ..., c\ are obtained 
in the vicinity of the threshold as: b\ = 0.03, b 2 = 0.6, C\ = —1.2, c 2 = 1. Notice that 
the amplitude C\ controlling the decay of the part of the wave function, which is due to the 
semiflexible nature of the polymer, differs from zero at the localization transition. 

Let us compare the adsorption of the semiflexible polymer with that of a flexible chain. 
The threshold value u t ji for adsorption of a flexible polymer onto a surface is obtained as 
4.94 in units used for the semiflexible polymer (w = 1), and is thus smaller than the value 
Ut = 7.8 obtained for the semiflexible polymer. We have obtained that both thresholds 
approach each other with increasing the width of the potential well w. The comparison of 
the monomer density as function of the distance to the surface for flexible and semiflexible 
polymer is plotted in Fig.2. Notice that the lengths and the energy are given in units stated 
in the text after Eq.(10). The distribution of the ends of the polymer, which is given by 
the wave function itself, behave qualitatively in the same way. Notice that the monomer 
density of the semiflexible polymer in the vicinity of the wall is lower than that for flexible 
polymer. This explains the higher value of the threshold for the semi-flexible polymer. To 
understand this qualitatively let us consider the effect of the wall on an ideal polymer coil, 
which is brought in a weak contact with the wall. It is intuitively clear that the deformation 
of the semiflexible coil demands higher external force, or equivalently the repulsion force of 
the wall is higher for semiflexible coil. The lower density of the monomers as a function 
of the distance to the wall for small distances is the consequence of this circumstance. 



Notice that our prediction is in disagreement with MD simulations j7|, which predict that 
the semiflexible polymer adsorbs easily. The adsorption state was studied in under the 
condition that one end of the polymer, which has finite number of segments, is fixed at the 
surface. The adsorption threshold obtained in this way corresponds to the case, when the 
loop length becomes comparable with the contour length of the polymer under investigation. 
At the localization transition threshold we obtained, u = ut, the loop length is infinite. The 
adsorbed state studied in M is expected to correspond to the regime of strong adsorption, 
where an opposite behavior as we predicted is expected (see for a related discussion of 
adsorption onto an interface). 

As in the case of adsorption onto a symmetric interface, the wave function becomes zero 
at — e = 1, which imposes a restriction on the applicability of the present theory, which 
thus is restricted to the regime of weak adsorption. The value — e = 1 corresponds to the 
depth of the potential well u sa ^ w = 15.5 for w — 1. The effect of the rodlike character of 
the localized polymer is expected to be strong for u > u saw . This is in agreement with the 



prediction of a liquid-crystalline phase made by Kuznetsov and Sung [10 . To describe the 



strong adsorption (u > u sa , w ), the higher derivatives in Eq.(|9]) have to be taken into account. 

V. CONCLUSION 

To summarize, we have considered adsorption of a semiflexible polymer chain onto inter- 
faces and surfaces by using the differential equation for the end-to-end distribution function 
of the polymer chain associated with the moment expansion of the latter. The present study 
shows that the adsorption of the semiflexible polymer is qualitatively different from that of 
a flexible one. The essential features of adsorption of the semiflexible polymer are: (i) the 
existence of a new local length scale, which results in two-exponential decay of the monomer 
density of adsorbed polymer; (ii) the binding energy of the semiflexible polymer is weaker 
than that for flexible polymer for both interface and wall. The approximative theory pre- 
sented here is restricted to the regime of weak adsorption (u < u sa>w ), where the effect of 
the rodlike behavior of the polymer on small scales is weak. 
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Figure captions 



Fig.l The interaction potential with the interface and the wall. 

Fig. 2 The monomer distribution ty(z) 2 of adsorbed flexible (dotted line) and semi-flexible 
polymer at a surface. Uq — 9, w — 1. 
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